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1 Introduction 

In the Standard Model (SM) the decays B® d -^ are mediated by flavour chang- 

ing neutral currents (FCNC). They are CKM suppressed, and happen via higher-order 
electroweak Feynman diagrams. Their decay rates, additionally, suffer helicity sup- 
pression. For the B® — > n + fi~ channel a rate of (3.54±0.30) x 10~ 9 has been predicted 
[D[2], and for B% ->- /i+ fi~ the prediction is (1.07 ± 0.10) x 10" 10 [TJ. In Standard 
Model extensions these rates may be enhanced. Any deviation from the SM pre- 
dictions will indicate new physics. These channels are important probes for physics 
beyond the SM. 

This paper describes the latest measurement^! of B® d — > from the ATLAS 

and CMS experiments [31 H] at the LHC. 

2 Analysis Strategy 

The branching fractions BR(B® d — > can be measured relative to a well mea- 

sured reference channel to minimize uncertainties, e.g. detector acceptance uncer- 
tainties, while keeping the analysis independent of luminosity variations and bb pro- 
duction cross-section uncertainties. For this purpose both ATLAS and CMS use 
B^ — > J/ipK — > [i + fi~K ± as the reference channel, because of large available statis- 
tics, and similar (di-muon) final states. To keep systematic uncertainties low, similar 
selection cuts are applied to both the signal and the reference channel events. 
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The branching fractions are computed as 



BR{B° sd -»■ fi + fi-) = / — x BR{B ± J/tfrK* -»■ n+^K^) x x (1) 

Nj/ipK± Js,d 

where N^+^- and Nj/^ K ± are the number of observed signal and reference channel 

events, respectively. The factor Rae = A -' /i ' K±£j/ ^ K± [ s to correct for the detector 

acceptances (A) and event selection efficiencies (e) estimated for the two channels 
using Monte Carlo events (MC). The ratios f u / f s ,d are the ratios of the 6-quark 
hadronization probabilities to correct for the different production rates of B and 
B° sd . The reference channel branching fraction BR(i? =t — y J/il)K ± — y y&jjTK^ 1 ) = 
(6.01 ± 0.21) x 10" 5 , and the ratio f s /f u = 0.267 ± 0.021 are taken from other 
measurements [51E]- The ratio fd/f u is taken to be 1 [TJ. 



3 ATLAS Analysis 

The ATLAS analysis expresses the branching fraction as a product of the observed 
number of signal events and a Single Event Sensitivity (SES): 

BR(B° S -»■ /I+/J,-) = N^- x SES. (2) 

For a single observed signal event, the branching fraction BR(5° — y would be 

equal to the SES. 

The analysis uses pp collision data at yfs = 7 TeV recorded by the ATLAS detector 
in the period April-August 2011. This corresponds to an integrated luminosity of 
2.4 fb _1 . The analysis is robust against pileup effects. The details of the analysis can 
be found in reference [8]. 

A topological trigger selects di-muon candidates above a transverse momentum 
(Pt) threshold of 4 GeV. The signal channel events in the invariant mass range 
m ^+ii- ^ [5066, 5666] MeV are hidden in the analysis ('blind' analysis) until event 
selection cuts are finalized. The events in the side bands {m^+^- G [4766, 5066] MeV, 
"V+/x- ^ [5666, 5966] MeV) are split into two sets. To avoid statistical biases, one 
set is used for cut optimization (odd numbered events in data), and the other for 
estimation of the background under the signal (even numbered events in data). 

The event selection cuts are optimized using a multivariate technique. The method 
uses Boosted Decision Trees (BDT) with 14 input variables to compute an event 
classifier, Q. The signal search window, Am A1 + M - , in the B° s —y fi + fi~ invariant mass 
spectrum is optimized together with the BDT classifier Q. The method determines 
the optimal Am M + M - and Q to get the maximum value for the estimator: 

V{Q,Am li+ll -) = "7==> (3) 

1 + J N bkg 
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where e S i g is the signal selection efficiency and Nf,k g is the continuum background 
interpolated from the side bands (two times the number of odd events). 

The events are split into three mass resolution categories distinguished by the 
maximum pseudorapidity (|^ M | ma a;) of the muon tracks. The Am (J + fl - and Q are 
separately optimized for each category. The same classifier cut is used to compute 
the acceptance and efficiency ratio R^ e from the MC, and the reference channel yield, 
Nj/^ K ±, in the three categories. The yield is determined by fitting the — >• 
J/%l)K ± — > ji + n~K ± invariant mass spectrum, and computing the Nj/^ K ± in a band 
mj^;c± e [5180, 5380] MeV. The SES is thus computed for the three categories using 
Equations (CQ) and (J2]). In the di-muon invariant mass spectrum there is an expected 
background due hadronic decays, where hadrons are misidentified as muons. The 
B —7- hh background is the dominating resonant background. It is estimated from 
the MC. 
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Figure 1: The di-muon invariant mass spectrum (points) (in the ATLAS analysis 
[8]) after unblinding in the three mass resolution categories. The search window is 
shown by dashed blue lines. The shaded grey regions reflect the sidebands used in 
the analysis. The solid curve shows the B° s — > fi + n~ distribution in MC scaled by a 
factor of 10. 
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3.1 Results 



The upper limit on the B s — > / u + pT branching fraction is computed using the CL S 
method [9]. A likelihood expression combines the SES computed in the three mass 
resolution categories, its uncertainties, and the expected resonant (B — > hh) and 
non-resonant background contributions to the B° s — > n + fi~ invariant mass spectrum 
[8]. An expected limit of (2.3^:°) x 10~ 8 on BR(B° S is obtained. The data 

in the blinded region is analysed, and the signal yield is measured in the optimized 
search window for the three mass resolution categories (see Figured]). The observed 
upper limit is 2.2 (1.9) x KT 8 at 95% (90%) confidence level (CL). The observed limit 
is comparable with the expected limit (Figure [2]). 
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Figure 2: The observed CL S (points) as a function of B° s — > branching fraction 

in the ATLAS analysis [8J. The upper limit is read at the intersection of the red 
line with the observed CL S line, which corresponds to 95% CL. The green and yellow 
bands indicate ±lcr and ±2cr deviation from the expected limit (dashed line). 



4 CMS Analysis 

In this section the key features in the CMS analysis different from the ATLAS one 
are presented. 

The CMS analysis uses 5 fb -1 of pp collisions at yj~s = 7 TeV recorded by the 
CMS detector in the year 2011. In this analysis a cut-and-count approach is taken. 
The optimization of the selection cuts is performed using the signal MC and all the 
side-band (m M+M - e [4900,5200] MeV, m M+M - G [5450,5900] MeV) events in the 
data. The events inside the signal region (in data) are kept blinded until the cuts are 
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established. A random-grid search method is used to tune the cuts on 11 analysis 
variables to get the best upper limit [11]. A different set of selection cuts is used for 
events in the endcap region of the detector than in the barrel region. 

The analysis shows good agreement between the reconstructed distributions in 
data and MC, and it is not sensitive to pileup events. 

4.1 Results 

The branching fractions for both B d — > fi + fi~ and B® — > fi + fi~ decays are measured 
simultaneously using two asymmetric search windows around the B d and B® masses. 
Figure [3] shows the di-muon invariant mass spectrum after the unblinding. 
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Figure 3: The di-muon invariant mass spectrum (in the CMS analysis [TT]) after 
unblinding in the two detector regions. The solid and the dashed horizontal lines 
indicates the search windows for B° s and B d events, respectively. 



The limit extraction takes into account the expected combinatorial and the res- 
onant background contributions (estimated using MC), and also the number of ex- 
pected signal events assuming the SM branching fractions. The resonant background 
includes B® d — > h+h^^ decays, where hadrons are misidentified as muons, and 
B® d — > h~n + v rare semileptonic decays. The limits are computed using the CL S 
method [HJC3II]- The expected upper limits for B° s — >■ (B d — > are 8.4 x 10~ 9 

(1.6 x 10" 9 ) at 95% CL [II]. The upper limit for B° d ->• n+fi' is 1.8 (1.4) x 10~ 9 at 95% 
(90%) CL. Figure 0] shows the dependence of CL s+ b on the B® — > branching 
fraction. 
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Figure 4: The observed CL s+b (solid line) as a function of B Q S — > ji + /i~ branching 
fraction (in the CMS analysis [11]) assuming the SM expectations. The green and 
yellow bands indicate ±1<t and ±2<r deviation from the expected limit (dashed line). 

5 Summary 

The B® — > fi + fi~ and B® — > are important decay channels to search for new 

physics. The latest limits on their branching fractions were set by the LHC ex- 
periments. These are summarized in Table [TJ The best B® — > fi + fi~ upper limit, 
4.2 x 10 -9 , is the combined result [13] of the measurements performed by ATLAS, 
CMS and LHCb experiments. It is close to the SM expectation, (3.54 ±0.30) x 10 -9 . 
The ATLAS and CMS collaborations are working on extending their analyses to the 
full available statistics collected from the LHC. 
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Table 1: Limits at 95% CL measured by different LHC experiments [81 fTT l IT2l [T3] . 
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